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00 

Let  f e L^(0,“>),  6 > 0 and  (G  f)  (t)  = 6 ^ / e^^  ^^'^^f(s)ds.  Given 

^ t 

a partition  P = {0  = t^  < t^  < • • • < t^  < < •••}  of  [0,<”)  where 

t^  ->•  oo,  we  approximate  f by  the  step  function  Apf  defined  by 

A f(t)  = (G.  G G.  f ) (0)  for  t,  < t < t.  , 

P 0 . 6 . , 6,  1-1  — 1 

1 1-1  1 

vrtiere  main  results  concern  several  properties  of  this 

■ process,  with  the  most  important  one  being  that  Apf  -»■  f in  L^(0,“>)  as 

u(P)  = sup  6.  ->  0.  An  application  to  difference  approximations  of  evolution 

• i 


SIGNIFICANCE  AND  EXPLANATION 

A method  of  approximating  functions  which  are  integrable  on 
by  piecewise  constant  functions  is  presented  and  studied 
in  this  paper.  The  method  used  and  the  properties  established 
for  it  allow  one  to  reduce  the  study  of  the  convergence  of  a 
difference  method  of  theoretical  interest  for  nonlinear  time- 
dependent  problems  with  forcing  terms  to  the  simpler  study  of 
related  problems  without  forcing  terms. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descrip- 
tive summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


AN  APPROXIMATION  OF  INTEGRABLE  FUNCTIONS  BY  STEP 


FUNCTIONS  WITH  AN  APPLICATION 


M.  G.  Crandall^^  and  A.  Pazy^^'^^ 


This  note  is  concerned  with  an  interesting  method  of  appr  ximating  an  integrable 
function  f : (0,®)  -*■  R by  step  functions.  The  approximation  pr^icess  involves  the 
integral  transformation  G^  : L^(0,“)  -*•  L^(0,”)  defined  for  6 > 0 by 


(1) 


(G^fXt)  = i/ 


e<'^-='/*f(s)ds 


Equivalently,  g *-  G^f  is  the  unique  function  g c L (0,»)  which  satisfies  g - 6g'  = f. 
Let 

P = {0  = t„  < t^  < ...  < t.  < t.^^  < ...) 

be  a partition  of  t0,“)  with  lim  t.  = “.  The  step  sizes  of  the  partition  are  denoted 

i-x»  ^ 

by  ~ Each  partition  P determines  a piecewise  constant  approxima- 

tion Apf  of  f defined  by 


(2)  Apf(t)  = (G,  G^ 

1 1-1 


G.  f)  (0)  for  t.  , < t < t.  , 
0^  1-1  — 1 


i = 1,2,...  . 


The  mesh  of  the  partition  is  denoted  by  p(P);  p (P)  = sup  6..  The  main  results  are 

l<^i<«)  ^ 

suranarized  in  the  following  theorem. 

Theorem:  Let  P be  as  above,  f e L^(0,“>)  and  A^  be  defined  by  (2).  Then 


(3) 


(4) 


(5) 


and 


(6) 


Apf  e L (0,®)  , 


J |a  f (s)  [ds  <_  / |f(s)lds  , 


J A f(s)ds  = / f(s)ds  , 


lim  / |a  f(s)  - f(s)|ds  = 0 
y(P)^-0  0 
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The  definition  of  the  transformation  f -►  A f as  well  as  the  questions  resolved 

P 

by  the  theorem  arose  naturally  from  considering  difference  approximations  of  certain 
nonlinear  evolution  problems.  While  this  motivation  is  not  relevant  for  the  statement 
or  the  proof  (given  in  Section  1)  of  the  theorem,  we  do  explain  it  briefly  in  Section  2. 
We  are  indebted  to  Carl  de  Boor  for  his  advice  on  this  problem. 


Section  1.  (Proof  of  the  Theorem) . 

Let  h,k  e L^(-",0)  and  f e L^(0,<»).  We  define  )c  » f e L^(0,“>)  and 
h * )c  e L^(-<«>,0)  according  to 


(k  » f)  (t)  = / k(t  - s)f(s)ds 


h » k(r)  = / h(r  - s)k(s)ds  . 


The  convolution  operator  is  commutative  and  associative,  while 


h o (k  o f)  = (h  * k)  « f . 


For  6 > 0 we  set 


The  transformation  G in  (1)  is 
0 


k^(r)  = 6 exp(r/6)  . 


Let  P={0=t„<t,  < •••<t.  <t.  , < •••}  and  6.  = t.  - t.  , be  as  in  the  introduc- 

01  1 1-1  1 1 1-1 

tion  and  be  given  by  (2) . For  simplicity  of  notation  we  will  set 

(1.6)  k.  = k,  and  K.  = k.  * K.  , = k.  * k.  , * •••  * k for  i = 1,2,...  . 

10  11  1-1  1 1-1  1 
i 

Since  k,  > 0,  A clearly  satisfies 
1 — P 

(1.7)  |Apf|  lAp|f|  . 

Moreover  by  (2),  (1.5),  (1.6)  and  (1.3) 
t. 

1 1 

/ Ap|f|(s)ds  = o (k^_^  o (...  „ (k^  „ (k^  „ |f  |))...)))(0) 


= I « |f|)  (0)  = / I 6 K (-S)  |f(s)  Ids 


i 

X 


( 

I 

J 


1 


since  each  of  the  summands  6j^Kj^(-s)  in  the  last  integrand  is  nonnegative,  we  can 


establish  (3)  and  (4)  of  the  theorem  by  showing  that 


(1.9) 

while  (5)  requires 


^ 1 for  -<»  < r < 0 


t=l 


(1.10) 


Y 6^K^(r)  =1  a.e.  on  -“  < r < 0 . 


li=l 


The  following  lemma  implies  (1.9)  since  » 1 ^ 0. 


Lemma  1 ; For  each  j = 1,2,... 

(1.11)  I + K.  * 1 = 1 


i=l 


Proof  of  Lemma  1:  We  proceed  by  Induction.  If  j = 1 the  claim  is  that  + )tj^  * 1 : 1. 


Indeed,  for  any  6, 

(1.12) 


r 


6)c^  + )t.  • 1 = e ' ” + T J e'‘  "'“ds  = 1 . 
0 0 0 


We  now  assume  the  claim  is  true  for  j = i and  verify  it  for  j = i + 1.  By  (1.12)  we 
have 


i+1 


(1.13)  K.  *1  = )t.  *K.*l  = )i.  •l.K.  = (l-6_  k )*K.  =K.»1  - (S.^  = 1 - I 6 K„ 

1+1  1+1  1 1+1  1 1+1  1+1  1 1 1+1  1+1  8 I 


where  the  last  equality  follows  from  the  induction  hypothesis.  Rearranging  (1.13) 
yields  (1.11)  with  j = i + 1 and  the  proof  is  complete.  • 

We  verify  (1.10)  indirectly.  Let 

-ot 


(1.14) 

If  o > 0 


f^(t)  = e 


(1.15) 


'=6^0)  <»^>  = i / 


(t-s)/6  -os 


e ds 


1 -ot 
e 


1 + o6 


from  which  it  follows  that 


(1.16) 


(ApfgXt)  = TT  *^i-l  - ‘ ^ *^1 


-3- 


i 


1 


and  hence 


(1.17) 


/ (Apf^)  (s)ds  = [ 6.  TT  (1  + 


i=l 


Setting  = <1  + a4^ 


-1 


we  claim  that 


(1.18) 


5 r,  + 6„r„r,  + ...  + 6.r.r,  , 
11  221  11  i-1 


i=l 


*^2^1  ^ ’^i’^i-1 


"2^1  = ° 


for  i = 1,2,...  . The  proof  parallels  the  proof  of  Lemma  1.  Since  ^ 6,  = ", 

i=l  ^ 

we  have  2.  ° i " and  (1.17),  (1.18)  together  imply 


/ (A„f  ) (s)ds  = i = / f (s)ds  . 

^ p 0 ft  ^ rr 


a a 


Setting  f = f^  in  (1.8),  letting  i .*•  <»  and  using  the  above  implies  (1.10). 

It  remains  to  verify  (6) . (We  remark  that  the  previous  results  did  not  require 
U(P)  < <»)  . In  view  of  (4),  which  is  independent  of  P,  it  suffices  to  verify  (6) 
for  a dense  subset  F of  L^(0,‘»).  It  is  convenient  to  choose  F = span{f^  : o > 0). 
(In  fact,  span{e  : n = 1,2,...}  is  dense  in  L^(0,'“),  as  is  well  known.  To  see 

this,  use  the  change  of  variables  x = e ^ which  exchanges  (0,")  and  (0,1)  while 
e becomes  x". ) To  proceed,  we  estimate  ~ terms  of  u(P).  For 

convenience  of  future  referencing  the  simple  lemma  which  does  so  is  stated  without  using 
the  notation  a>x>ve. 

00 

Lemma  2;  Let  ® sequence  of  positive  numbers  satisfying  ^ 6.  = “ and 

i=l 


o > 0.  Let  t„  = 0,  t.  = 6,  + 6,  + 
0 1 1 2 


+ 6.  for  i = 1,2,...  and  y = sup  S..  If 

l<i<«j 


g(t)  = 


i"l 


(1  + aS^) 


for  t.  , < t < t, 
1-1  — 1 


then 


2 2 

I -oti  -ot  , uat  ay  , , -oy, 

lg(t)  - e I < e max{e  e - 1,  1 - e > . 


[ 


Proof  of  Lemma  2;  It  is  enough  to  treat  o = 1,  for  then  the  general  result  follows 
upon  replacing  {6^^}  by  {06^)  and  t by  ot.  Elementary  calculus  yields 

,2 


(1.19) 


e < (1  + 0)  < e 


for  6 > 0 


-4- 


Multiplying  these  inequalities  yields 


-(6  +...+6.)  i , -(6, 

2 "■  < TT  (1  + 1 e 

£=1 


I e 


Using  (1.20)  and  the  inequalities 


6^  + •••  + 6?  < p(6  + •••  + S.)  and  t<6,  +-*'+6.  <t+u 

1 1—1  1 — 1 1 — 


for  t.  , < t < t.  establishes 
1-1  — 1 


- 1)  < g(t)  - e'^  < - 1) 


and  hence  Lemma  2 in  the  case  0=1. 


End  of  proof  of  the  Theorem:  By  (1.16)  and  Lemma  2 


/ [Aj-f  - f Ids  < J e °^max{e*^*^^e'^  - 1,  1 - e °^}ds 

0 P ° ° ~ 0 

where  y = y(P).  The  right  hand  side  above  tends  to  zero  as  y -♦  0,  and  the  proof  is 
complete.  • 

1 oo  -1 

Remar )t : If  m £ L (-",0)  n l (-“,0),  = 6 m(r/6)  and  (1)  is  replaced  by 

Ggf  = m^  o f,  the  first  part  of  the  above  proof  adapts  to  establish  that 


/ |A  f(s)  |ds  < C / |f  (s)  |ds 


provided  C ^ 0 can  be  chosen  to  satisfy 


i(r)  I + C / |m(s)  |ds  ^ C a.e.  -<»  < r £ 0 . 


This  last  estimate  is  equivalent  to  |m.  | + |m,  [ * C ^ C.  If  ||m|| 

L^  (-«’,0) 


< 1 we  may  set 


I < 


/(I  - m 


L (-»>,0) 


L (-”,0) 


(Consideration  of  the  case  m(r)  = 2e  shows  some  such  restriction  is  necessary.)  By 


contrast,  our  proofs  of  (5)  and  (6)  required  special  properties  of  the  exponential  )cernel. 


U'  + AU  3 0 

(2.4)  ] 

[ y(0)  = (x,f)  . 

By  the  (clearly  valid)  version  of  Theorem  1 in  which  L^tO,™)  is  replaced  by 
the  term 

N-1  *^i+l 

I / <0)  ■ f (s)  ll'^s 

i=0  t. 

1 


tends  to  zero  as  u = max  6.  -►0. 

1 

i 

We  now  briefly  sketch  the  proof  that  A is  accretive.  If  Z is  any  Banach  space,  set 


Up  + ^<1112  - IIPII2  Up  + ■ IIpIU 

Ip<ql„  = = inf  

^ AlO  A A>0  A 

If  Y = X X L^(0,";X)  is  equipped  with  the  norm 


(x,g)  lly  = l|x||^  + / ||g(s)||jjds  for  (x,g)  e x 


then  one  computes  that 


I (x,g) , (y,h) ] = Ix,y)  + / Ig  (s) ,h  (s) 1 ds  . 

y X ^ X 

It  follows  that  A is  accretive  in  Y if  for  every  (x,y) , (x,y)  e A and 
g,g  e W^'^(0,";X)  we  have 

(2.5)  [ (x  - x,g  - g)  , ( (y  - g (0)  ) - (y  - g (0)  ) , -g'  + g’ ) 1 ^ 

00 

= [x-x,(y-y)  - (g(0)-g(0))l  + / [g(s)-g(s),-g'(s)+g'(s)  1 ds  > 0 . 

The  fiist  term  in  (2.5)  is  estimated  by 

(2.6)  [x  - X,  (y  - y)  - (g(0)  - g(0) ) I > lx  - x,  (y  - y)  - ||q(n)  - g(0)  |1  ' !lg(P>  ' g(0>  II. 

X X A 

where  the  first  inequality  is  due  to  the  fact  that  [p.ql^  is  Lipschiiz  in  q with 
constant  1 and  the  second  inequality  is  because  (x,yl,(x,v)  r A and  A is  accretive. 
The  second  term  in  (2.5)  is  given  by 


(2.7)  / [g(s)-g(s),-g’(s)*g’(s)l  ds  = / - — ||g(s)-g{s)|l  ds  = ||q(0)  - q(n)|L 

0 X Q ds  X X 


-7- 


since  ^ l|k(s)lljj  = Ik (s)  ,k'  (s)  1 wherever  l|k(s)||jj  and  k(s)  are  both  differentiable. 
Together,  (2.6)  and  (2.7)  imply  (2.5)  and  hence  that  is  accretive. 

The  system  (2.4)  was  introduced  in  [4]  for  another  purpose.  The  results  we  have 
just  obtained  reduce  many  questions  concerning  (2.1)  to  the  same  questions  for  f = 0. 
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